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Abstract. We study the possibility of consistent extension of MDI full-disc helioseismic 
campaigns with the growing data set of HMI observations. To do so, we down-sample and 
filter the HMI Dopplegrams so that the resulting spatial power spectrum is similar to the 
spatial power spectrum of MDI full-disc Dopplergrams. The set of co-spatial and co-temporal 
datacube pairs from both instruments containing no missing and no bad frames were processed 
using the same codes and inverted independently for all three components of the plasma flow 
in the near surface layers. The results from the two instruments are highly correlated, however 
systematically larger (by ~ 20%) flow magnitudes are derived from HMI. We comment that this 
may be an effect of the different formation depth of the Doppler signal from the two instruments. 



1. Motivations 

Understanding the details of the solar dynamo is inevitably linked to measurements of plasma 
flows inside the Sun. Methods of helioseismology [U [2] can measure such flows, one of them 
being time-distance local helioseismology [3j. It is based on the measurements of travel times 
of seismic waves, which propagate through the solar interior and carry information about the 
inhomogeneities (flows in particular) affecting their travel time. By solving the inverse problem 
one can learn about these inhomogeneities. 

A proper assessment of the character and structure of solar flows requires a long term study. 
Recently, the Michelson Doppler Imager (MDI [5J, on board of SOHO spacecraft), which covered 
15 years of observations, was replaced by its successor Helioseismic and Magnetic Imager (HMI 
[5], on board of SDO spacecraft). Our study is thus driven by the need to fully exploit all of 
the data. Here, we specifically ask the question: Can MDI and HMI data be possibly utilised 
together for inverse modelling for time-distance local helioseismology? 

In local helioseismology, two methods are usually used to solve the inverse problem. The 
regularised least squares (RLS; in time-distance helioseismology used for the first time by [6]) 
minimisation seeks to find models of solar interior that minimise the difference between the 
predicted and measured travel times. The inverse problem must be solved separately for each 



Figure 1. Cuts through averaging kernels 
of inversion for horizontal v x and vertical v z 
flow in the middle point as a function of 
depth. Note that z = at an optical depth 
of T500 = 1, negative z values refer to solar 
interior and positive ones refer to atmosphere. 

set of travel times. On the contrary, optimally localised averaging (OLA; 0IHJI9]) method is 
based on explicitly constructed spatially confined averaging kernels, while bounding the noise 
propagation at the same time. The inverse problem in the OLA approach must be solved only 
once and does not directly involve the measurements, except for the knowledge of statistical 
properties of the travel-time noise. Estimates for inverted quantities are then obtained by 
taking linear combination of travel times utilising a set of weights resulting from the inversion. 
When the measurements from various instruments are combined to form a homogeneous data 
set, the OLA is a natural approach to be used for inverse modelling from such data set. 

To test inverse modelling, we constructed a well-understood inversion for all three flow 
components. The inversion follows the Fourier-space based solution [10] using the validated 
inversion code implemented within the German Science Center for SDO at Max-Planck- 
Institut fur Sonnensystemforschung in Katlenburg-Lindau, Germany. This code introduces 
additional terms in the cost function of the inverse problem, which allows the minimisation of the 
contributions from cross-talks. This property is crucial in order to be able independently invert 
also for a weak vertical flow. The inversion utilised a combination of wave sensitivity kernels 
computed in the Born approximation |12| for the surface gravity wave (/) mode, assuming the 
application to the MDI full-disc Dopplergrams. The travel-time noise covariance matrix used in 
the inverse problem was measured from the data using 1/T fitting approach |13j . The inversion 
combined sensitivity kernels for travel times for three point-to-annulus and point-to-quadrant 
geometries (making the measurements sensitive to divergence of the flow and vertical velocity, 
and to waves travelling in west-east and north-south directions) with radii of the annuli from 
7.3 to 29.2 Mm. The inversion targets only near sub-surface layers (see the averaging kernels in 
Fig.0. 

2. Data 

To do a proper comparison between the two instruments, we required the consistent analysis of 
the data, which are both co-spatial and co-temporal. MDI delivered data suitable for local 
helioseismology only two months a year (during Dynamics Campaigns), with the last one 
(June and July 2011) overlapping with routine operations of HMI. In this two- month period 
we prepared 64 consecutive pairs (one for MDI, one for HMI) of 12-hour long Dopplergram 
datacubes capturing the centre of the solar disc using the mapping and tracking pipeline at the 
GDC-SDO. Additionally, we did not want to have results affected by missing or bad frames in 
datacubes, thus we required a 100% duty cycle. A down-selection was done on a frame-by-frame 
basis manually. Only 4 (!) out of 64 co-spatial and co-temporal datacube pairs fulfilled our 
requirements. 
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Figure 2. Left: Spatial power spectra of MDIfd Dopplergrams compared to HMI Dopplergrams 
down-sampled to MDIfd resolution before and after correction. Right: The wave-number filter 
designed for the correction. 

The sensitivity kernels used in the inverse problem described above were designed to be 
consistent with MDI full-disc (MDIfd hereafter) Dopplergrams. HMI is a higher-resolution 
instrument, thus we needed to make sure that HMI data looked like MDIfd ones. This was 
achieved by a Postel-projection of both MDI and HMI Dopplergrams with the same map scale 
(corresponding to MDIfd pixel size) using bicubic interpolation. HMI Dopplergrams were thus 
down-sampled to MDIfd resolution. Even in such data the contribution from large wave-number 
was much higher than in the case of true MDIfd data (see Fig. [2] left) and we wanted to correct 
for it. This may be done by an application of a filter in the wave-number space (Fig. fright), 
designed so that the spatial power spectra of both instruments decay similarly in the large 
wave-number region. This filter approximates the correction to the different modulation transfer 
functions (MTF) of the two instruments. Note that despite the filtering, spatial power spectrum 
of corrected HMI Dopplergrams still have an excess in power at angular degrees I of 50 to 300. 
We will return to this issue in the next section. 

Apart from the power-spectrum correction, the rest of the analysis (filtering, travel-time 
measurements, application of inversion weights) was identical for both MDI and degraded HMI 
datacubes. 

3. Comparison of helioseismic observables 

Consistently with the setup of the inversion, we measured wave travel times from /-mode- 
pass filtered Dopplergram datacubes with point-to-annulus and point-to-quadrants geometries 
with distances 7.3 to 29.2 Mm. Travel-time maps measured from corresponding MDI and HMI 
datacubes are highly correlated (correlation coefficients of 0.8 or higher), but the travel times 
measured from HMI datacubes are systematically larger by 20% (e.g., Fig.|3]left). This statement 
holds for all pairs of travel-time maps. The excess is mostly prominent on supergranular and 
similar scales (angular degree between 50 and 300) and it seems to be a consequence of the 
excess power discussed above (see Fig. [2] left). 

Since the estimate of the inverted flows is just a linear combination of travel times and the 
inversion weights are the same for both MDI and HMI measurements, the magnitudes of flow 
estimates are also systematically larger, when the flow is inverted from the HMI datacubes than 
from the MDI datacubes (see example in Fig. fright). Apart from this systematic effect (there 
is an excess of 20% for horizontal components and around 30% for vertical component) the flows 
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Figure 3. Left: Scatter plot of the travel times measured for the /-mode with point-to-annulus 
geometry with distance 10 px from MDIfy (horizontal axis) and degraded HMI (vertical axis). 
Right: Equivalent scatter plot to compare the inversion estimate of v x with the correlation 
coefficient indicated. In both cases, the red solid line represents the best least-squares fit of the 
points and the black solid line has a unity slope to increase the visibility of the systematically 
larger values derived from HMI Dopplergrams. 

derived from MDI and HMI measurements are almost identical (correlation coefficient larger 
than 0.9 for horizontal components and around 0.8 for the vertical component). The lower 
correlation in the case of vertical component is probably caused by lower signal-to-noise ratio of 
the vertical flow estimate. 

An example of full-vector flow maps for HMI and MDI is displayed in Fig. |4| One can 
see that except for the systematically larger magnitudes in the HMI flow map, both maps are 
practically identical. In these maps, supergranular cells (rosette-like features in the horizontal 
components) are dominant. Also note that the up flows (red) in the vertical component in most 
cases correspond to outflows in the horizontal component. The correlation between the vertical 
component and the horizontal divergence of the horizontal flow is higher than 0.8. This is one 
indication that the inversion for the vertical flow is reasonable. 

4. Final remarks 

We showed that with some drawbacks, HMI and MDI helioseismic datasets may be analysed 
together. Long-term studies of flows in the solar convection zone employing combined datasets 
are therefore possible. However, the noticed bias must be understood and corrected for. 

It is likely that the different flow magnitudes are an effect of the different spectral line used in 
the two instruments. The formation depth of HMI Doppler signal measured in Fel (617.3 nm) 
line is some 25 km deeper in the photosphere than the formation depth of the MDI Doppler 
signal measured in Nil (676.8 mm) line [Hj. Thus the convection (especially supergranules) is 
better pronounced when observed in the Fel line, hence the power spectrum is expected to have 
more power on scales where convection is dominant, including the supergranular scales around 
angular degree of 120. That is the case and it can be clearly seen in Fig. [2} Increased sensitivity 
to the convection affects the sensitivity kernels, which slightly differ for the two spectral lines 
(however, we did not compute the proper sensitivity kernels for Fel line, so we did not model 
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Figure 4. Comparison of full-vector flow maps in the same region on the Sun, inverted from 
MDI (left) and HMI (right) Dopplergrams. Horizontal components in the plane of solar surface 
are indicated by arrows, the vertical component by colours. The level of random error in the 
flow estimates is 49 m/s for horizontal flow components and 4 m/s for the vertical one. 



the effect properly). As a consequence, we expect that when the weights from the inversion 
performed with Nil-line sensitivity kernels are applied to travel-time maps measured from Fel 
line, the true inversion averaging kernel is slightly different, with the peak sensitivity located 
probably a little deeper (the averaging kernels for the Nil-line inversion are plotted in Fig. [TJ. 
If we assume that the magnitude of the plasma flow increases with height in the last 1 Mm of 
the convection zone and then drops gradually around T500 = 1 level, the averaging kernel with 
a deeper peak averages the flow with slightly higher amplitudes and therefore the resulting flow 
estimate will also have larger magnitude. That is observed as a bias when comparing the flow 
estimates from the two instruments. The higher bias in the vertical flow (~ 30%) than in the 
horizontal components (~ 20%) might indicate a steeper gradient of the vertical flow below the 
surface. 

Other effects may also be responsible for the excess power in the spectrum, such as a more 
complicated instrumental MTF shape which was not corrected for. After the bias is understood, 
it would be possible to correct for it, for example by application of an additional non-monotonic 
wave-number filter. 

This study is the first to compare inverse modelling using the OLA technique for MDI and 
HMI data. Similar direct comparisons of flow and sound speed inversions was performed on 
early HMI data (J Zhao, private communication) using the RLS technique, but it showed larger 
discrepancies not only in magnitude, but also in structure of the flows than presented here. This 
higlights the necessity to first test and understand the systematics due to both the instrument 
and the analysis technique. 
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